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======================================================================

** Conservation of Energy (system view)
*** System kinetic energy

K_sys = sum K_i = sum (1/2)m_i v_i^2

Partition theorem:
  K_sys = K_cm + K_rel
where
  K_cm = (1/2)M v_cm^2   KE associated with cm motion
  K_rel = sum_i (1/2) m_i v_ic^2  KE relative to cm
        = sum_i (1/2) m_i |v_i - v_cm|^2

For two-particle systems:
  K_rel = (1/2)mu v_12^2  (mu = (1/m1 + 1/m2)^-1 = reduced mass)
This special case is proved in an earlier video
  ["elastic parameter" in cons-p.txt]

Derivation:
  K_sys = sum_i (1/2) m_i v_i^2
        = sum_i (1/2) m_i |v_cm + v_ic|^2
        = sum_i (1/2) m_i (v_cm + v_ic) dot (v_cm + v_ic)
        = sum_i (1/2) m_i [v_cm^2 + v_ic^2 + 2 v_cm dot v_ic]
        = sum_i [(1/2) m_i v_cm^2]  +  sum_i [(1/2) m_i v_ic^2]
            + sum_i [m_i v_ic dot v_cm]
        = (1/2) [sum_i m_i] v_cm^2  +  sum_i [(1/2) m_i v_ic^2]
            + [sum_i m_i v_ic] dot v_cm
        = K_cm + K_rel + (0 dot v_cm)
        = K_cm + K_rel
[See "Momentum of a system" in cons-p.txt for derivation of P_rel = sum_i m_i v_ic = 0]

Note that K_rel is not necessarily zero -- KE can hide (contrast with
momentum, which cannot).
  [show page from "momentum of a system" in cons-p.txt]

K_rel can be further divided into rotational and vibrational
contributions.  More about that later.

======================================================================

*** Work done on systems
**** Work done by a force

The work done by a force on a system is basically the work done by
that force, never mind that it is part of a system:

[draw particles in a system -- with particle i acted on by F_i]

  W_i = int F_i dot dr_i

Note that the work done by a force is calculated based on the
displacement of the part of the system upon which it acts -- this
is VERY important.  Other parts of the system (including cm) may move
differently, so it is important to use the correct displacement.

W-KE theorem for single particle:
  W_i,net = dK_i

W-KE theorem for entire system:
  W_sys = dK_sys
where W_sys = sum_i W_i,net is the sum of the work done by every force
acting on every particle in the system.

As we will see later, it is important to include internal and
external forces.

======================================================================

**** Example [2-object spring mass from sys.txt]

We introduced this example when discussing systems of particles:

[see "Two masses connected by spring..." in sys.txt]
[see "Example: spring-mass oscillations with two masses" in sys.txt"]  <-- maybe not

Two blocks (m1 and m2) are connected by a massless spring (k) whose
natural length is L0 (dist between m1 and m2 centers).  Both blocks
are free to slide (no friction) across a horizontal floor.

At t=0, m1 is pressed against a wall (x1=0) and the spring is
compressed by d (x2=L0-d).  Both are released from rest at t=0.


Previously, we solved for the motion of the two blocks.  At this
point, we are interested in the work done by each force acting on
the system during the following time interval:

i = time t=0, when m2 is first released from rest
f = time t=t0, when m2 reaches equilibrium (x2=L0), at which
    point m1 will leave the wall

Reminder that, during this time interval:
  v1(t) = 0
  v2(t) = d*omega*sin(omega*t)
  vcm(t) = (m1*v1 + m2*v2)/M
         = (m2/M)*d*omega*sin(omega*t)
where
  omega = sqrt(k/m2)
  t0 = pi/(2*omega)

[Force diagrams:
m1: m1g, Fs12, N1G, N1W
m2: m2g, Fs21, N2G]

Four of these forces are vertical, hence perpendicular to any motion
that might exist:
  W_m1g = W_m2g = W_N1G = W_N2G = 0

Two of the remaining forces act on m1, which does not move during [0,t0]:
  W_N1W = W_Fs12 = 0

The remaining force is the spring force on m2, which we know how to
integrate:
  W_Fs21 = int Fs21 dot dr2 = (1/2) kd^2

The net work on the system is the sum of these results:
  W_sys = sum_i W_i = (1/2) kd^2

We can check the W-KE theorem by calculating K_sys,i and K_sys,f
  K_sys,i = 0  [neither block is moving]
  K_sys,f = (1/2)m1 v1f^2 + (1/2 m2 v2f^2
          = 0 + (1/2) m2 (d*omega)^2
          = (1/2) m2 d^2 (k/m2)
          = (1/2) kd^2
  DK_sys = (1/2) kd^2 = W_sys  [check]

It is interesting to note that the one force that does non-zero work
on the system is an internal force, and yet center-of-mass motion
changes, which implies a change in K_cm.  We will discuss this
shortly.

======================================================================

**** Center-of-mass and relative work

W-KE theorem for systems:
  W_sys = sum_i W_i = dK_sys  [work summed over all forces on all particles]

The work done by each force acting on a system can be divided into
its contribution to changes in K_cm and K_rel:
  W_i = W_i,cm + W_i,rel
  W_cm = sum_i W_i,cm = dK_cm
  W_rel = sum_i W_i,rel = dK_rel

  W_i,cm = int F_i dot dr_cm   [work done, pretending particle moves with cm]
  W_i,rel = int F_i dot dr_ic  [work done, based on displacement rel to cm]

Proof:

1.
  W_i,cm + W_i,rel = int F_i dot dr_cm + int F_i dot dr_ic
     = int F_i dot [dr_cm + dr_ic]
     = int F_i dot dr_i
     = W_i

2.
[see "Proof of Work-KE theorem" in cons-E.txt]
  dK_cm/dt = (d/dt) (1/2) M v_cm dot v_cm
           = M a_cm dot v_cm
           = sum_i F_i dot v_cm
Then
  DK_cm = K_cm,f - K_cm,i
        = int(ti->tf) (dK_cm/dt) dt
        = int(ti->tf) sum_i F_i dot v_cm dt
        = sum_i int F_i dot dr_cm
        = sum_i W_i,cm
        = W_cm

3.
  W_rel = W_sys - W_cm
        = dK_sys - dK_cm
        = dK_rel

======================================================================

**** Work done by internal forces (force pairs)

Work done by force on a system in general:
  W_i = int F_i dot dr_i

Forces can be divided into external and internal forces.  Same
formula applies to each individual forces, but for internal forces,
both forces in the force pair act on the system.  The combined work
done by the force pair is given by:

  W^int_ij = W_ij + W_ji
           = int F_ij dot dr_i  +  int F_ji dot dr_j
           = int F_ij dot dr_i  -  int F_ij dot dr_j
           = int F_ij dot (dr_i - dr_j)
           = int F_ij dot dr_ij   [dr_ij = displacement of i relative to j]

If we divide the work done by the force pair into cm and rel parts,
we find that
  W^int_ij,cm = int F_ij dot dr_cm  +  int F_ji dot dr_cm
              = int (F_ij + F_ji) dot dr_cm
              = int (0) dot dr_cm
              = 0

Internal forces do not contribute to W_cm.  Another way to see
this is to sum cm work due to ALL forces.  Since the displacement
is always dr_cm, the forces can be factored out and summed
individually -- the internal forces will cancel in pairs:

  W_cm = sum_i int F_i dot dr_cm
       = int (sum_i F_i) dot dr_cm
       = int F_ext dot dr_cm    [F_int cancels in pairs]

All of this is consistent with the idea that internal forces cannot
influence cm motion.


To summarize, we can divide forces acting on a system into two
categories: external and (always paired) internal.  Likewise, the
work done by each force can be divided into cm and rel contributions.

              cm          rel         total
external   W_ext,cm     W_ext,rel     W_ext  <-- transfer of energy in/out of system
internal      0         W_int,rel     W_int  <-- transfer of energy within system
all         W_cm         W_rel        W_sys
           = dK_cm      = dK_rel     = dK_sys

Note that W_int can increase/decrease K_rel (and therefore K_sys),
but that energy cannot come from the outside, since the forces doing
it are internal -- instead other forms of energy are involved.

======================================================================

**** Example revisited

Reminders:

Two blocks (m1 and m2) are connected by a massless spring (k) whose
natural length is L0 (dist between m1 and m2 centers).  Both blocks
are free to slide (no friction) across a horizontal floor.

At t=0, m1 is pressed against a wall (x1=0) and the spring is
compressed by d (x2=L0-d).  Both are released from rest at t=0.

We calculated the work done by each force over the following time
interval:

i = time t=0, when m2 is first released from rest
f = time t=t0, when m2 reaches equilibrium (x2=L0), at which
    point m1 will leave the wall

[Force diagrams:
m1: m1g, Fs12, N1G, N1W
m2: m2g, Fs21, N2G]

The net work is
  W_net = sum_i W_i = W_Fs21 = (1/2) kd^2  [all other W_i = 0]

This checked out with W-KE theorem.

Now let's calculate W_i,cm and W_i,rel for each force:

For m1g, N1G, m2g, and N2G, these forces are perpendicular to all
motion (individual particle, center of mass, relative), and thus
all work contributions are zero:
  W_m1g,cm = int m1g dot dr_cm = 0
  W_m2g,cm = W_N1G,cm = W_N2G,cm = 0  (same reason)
  W_m1g,rel = int m1g dot dr_1c = 0
  W_m2g,rel = int m2g dot dr_2c = 0
  W_N1G,rel = W_N2G,rel = 0  (same reason)

For Fs12 and Fs21 (internal forces, combined calculation), we find that
  W_Fs,cm = W_Fs12,cm + W_Fs21,cm = 0  (true for any internal force pair)
  W_Fs,rel = int Fs21 dot dr_21 = (1/2) kd^2  (dr_21 = dr_2, since m1 never moves)

N1W is interesting because it is an external force, does (net) zero
work, and yet W_N1W,cm and W_N1W,rel are both non-zero:

  W_N1W,cm = int N1W dot dr_cm
           = int (-Fs12) dot (m2/M)dr_2  [N1W=-Fs12 by N2L]
           = (m2/M) int (+Fs21) dot dr_2
           = (m2/M) (1/2) kd^2

  W_N1W,rel = int N1W dot dr_1c
            = int N1W dot (-dr_cm)   [dr_1c = dr_1 - dr_cm = 0 - dr_cm]
            = -W_N1W,cm
            = -(m2/M) (1/2) kd^2

Summary:

Force          W_cm                W_rel           total
----------------------------------------------------------
m1g (ext)        0                   0               0
m2g (ext)        0                   0               0
N1G (ext)        0                   0               0
N2G (ext)        0                   0               0
N1W (ext)  (m2/M)(1/2)kd^2   -(m2/M)(1/2)kd^2        0
Fs  (int)        0               (1/2)kd^2       (1/2)kd^2
----------------------------------------------------------
total      (m2/M)(1/2)kd^2    (m1/M)(1/2)kd^2    (1/2)kd^2

Let's compare W_cm and W_rel totals to K_cm,f and K_rel,f:

  v2f = d*omega   (omega = sqrt(k/m2))
  vcmf = (m2/M)d*omega
  K_cm,f = (1/2) M vcmf^2
         = (1/2) M [(m2/M)d*sqrt(k/m2)]^2
         = (1/2) M m2^2 / M^2 * d^2 * k/m2
         = (1/2) (m2/M) kd^2   [checks out]

  v21f = v2f = d*omega
  K_rel,f = (1/2) mu v21f^2   [mu = (1/m1 + 1/m2)^-1 = m1m2/M]
          = (1/2) m1m2/M d^2 (k/m2)
          = (1/2) m1/M kd^2  [checks out]

Note that not all kinetic energy ends up in the relative motion of the
two objects -- some of it gets transformed to K_cm (this occurs thanks to
N1W; we will discuss this mechanism in greater detail in the next video).

One consequence of this is that, once m1 leaves the wall, K_cm can no
longer change, and the amplitude of the spring-mass oscillations (in
terms of relative motion) gets reduced.

Recall that for t>t0, the relative motion is given by
  x21'(t) = x21(t) - L0 = A sin(omega'(t-t0))
  v21(t) = A*omega' cos(omega'(t-t0))   [show page from previous lecture]
  omega' = sqrt(k/mu)
  A = d*omega/omega'
    = d*sqrt(mu/m2)
    = d*sqrt(m1/M)

This results in a maximum relative velocity of
  v21_max = A*omega'   [occurs whenever U_spring=0]
which leads to
  K_rel,max = (1/2) mu (A*omega')^2
            = (1/2) mu A^2 (k/mu)
            = (1/2) kA^2
            = (1/2) kd^2 (m1/M)
which is consistent with K_rel at t=t0.

We will have more to say about this example when we cover system
mechanical energy.

======================================================================

**** External constraint forces

External constraint forces often do zero work on the system.  However,
there are circumstances where such a force can still play a role in
how energy is distributed between center of mass and relative motion.

The external forces which tend to do this are the ones that are
directed along the same line as the motion of the center of mass,
although they prevent motion of the part of the system upon which they
act.  In fact, the center of mass motion may very well be CAUSED by
this constraint.

Three examples...


Example 1: the two mass + spring example we just revisited

[force diagram]

There are three external constraint forces which do zero work in this
example: N1G, N2G, and N1W.

N1G and N2G do zero work because they act perpendicular to all motion
within the system (W_total = W_cm = W_rel = 0) -- we are not concerned
with these forces here.

The force we will discuss is N1W.  It does zero work on the system because
it acts on m1, which does not move during the time that N1W acts.

However, during the time N1W acts, m2 DOES move (thanks to the
spring), so the cm moves and W_N1W,cm is non-zero.

Examining the table of work values, one can see the impact of this
force.

Force      W_cm    W_rel    total
----------------------------------
m1g, etc    0        0        0

Fs          0      ++++++   ++++++   [actual energy source]

N1W         ++       --       0
----------------------------------
total       ++      ++++    ++++++

Without the wall, the system would be isolated, and the center of mass
would never move -- the spring would release its potential energy into
m1 and m2, and we would have a fair amount of relative motion, but no
cm motion.

With the wall in place, some of that energy is being transferred from
relative motion within the system to overall motion of the center of
mass.

This same principle applies in the N3L demostration that I showed earlier
[demo 050.m4v].  The muscles within my body do relative work on the system
containing me and the skateboard when I extend my hand.  The wall prevents
my hand from moving, and transfers some of the K_rel into K_cm, allowing
me to get moving.  Here again, the external constraint force in question is
a normal force from the wall.


Example 2: Car acceleration -- wheels

[see "Static friction between tires and road..." in forces.txt]

The second example involves a car accelerating from rest.  We discussed
this in an earlier video.

The car needs to be thought of as a system, since the wheels undergo
rotational motion, and therefore individual particles of the wheel
move relative to the cm of the car.

There is also relative motion within the engine of the car (pistons going
up and down, etc), which play a key role in providing an internal source
of energy needed to drive the car forward.

[force diagram: mg, N, fs(fwd)]

While the internal forces within the engine are responsible for
providing the KE to the car system, the force actually responsible for
accelerating the car is static friction between the tires and the
road.

Without friction, the tires would spin in place and the car would
not move -- plenty of K_rel, but no K_cm.

This static friction is an external force that does zero work:
  W_fs = int(ti->tf) fs dot v_contact dt
       = int(ti->tf) fs dot (0) dt
       = 0
[note: v_contact is the velocity of the part of the system in contact
with the road at a given instant -- this is the correct way to calculate
work done by fs when that force acts on different parts of the system
at different times]

However, fs is directed in the same direction as the motion of the
car's center of mass, and so
  W_fs,cm = int(ti->tf) fs dot v_cm dt  >  0
Likewise, the contact point on the tires are moving backward relative
to the axle (and hence the cm of the car), so
  W_fs,rel = int(ti->tf) fs dot v_contact,cm dt  <  0

The net effect of this static friction force is to transfer some
(actually most) of the relative work being done by internal forces
within the engine of the car (which facilitates wheel rotation) into
center-of-mass work done on the car as a whole.

The energy transfer is as follows:
  internal energy of gasoline / battery (electric cars)
     -> relative KE within the engine  [thermodynamic process]
     -> rotational KE within the wheels  [mechanical process]
     -> cm KE of the car  [static friction]

Force              W_cm     W_rel    total
-------------------------------------------
mg, N               0         0        0

internal forces     0       ++++++   ++++++  [actual source of energy]
within engine,
gears, etc

fs(tire/ground)    ++++      ----      0
-------------------------------------------
total              ++++       ++     ++++++


Example 3: disk rolling down an incline (static friction)

The third example involves a wheel rolling down an incline with
the help of gravity.  We will discuss this situation in detail when
we get to rotational motion later in the video series.

[demo 350]
[wheel on incline -- force diagram: N, mg, fs]

In the previous examples, the actual source of energy is internal to
the system, and is converted to KE by internal forces which do net
positive work on the system.  The external constraint force then
converts some of the resulting K_rel into K_cm.  In such cases,
W_cm>0 and W_rel<0 (while W_overall=0) for the constraint force
in question.

In the example we are considering now, the source of energy is an
external force (mg) that would normally produce a certain amount of
K_cm.

However, an external constraint force (fs in this case) will convert
some of the K_cm into K_rel, which in this context, is rotational KE.

Indeed, this static friction force is needed to ensure that the
rotational motion and translational motion of the wheel remain
properly corrolated so that the contact point is not slipping.

Indeed, v_contact=0 throughout this motion, and so W_fs overall is
zero (same calculation as in example 2)

In this case, W_fs,cm < 0 (fs and v_cm are in opposite directions) and
W_fs,rel > 0.  The static friction force is responsible for converting
some of the translational KE that is normally produced by gravity into
rotational KE of the wheel as the wheel accelerates down the incline.

We will have more to say about this situation when we cover rotational
motion.

Force    W_cm     W_rel    total
---------------------------------
 mg     ++++++      0      ++++++

 N        0         0        0

 fs       --       ++        0
---------------------------------
total    ++++      ++      ++++++

======================================================================

**** W-KE theorem in a non-inertial reference frame

W-KE theorems, including the breakdown into cm and rel contributions:
  K_sys = K_cm + K_rel;   W_sys = W_cm + W_rel
  W_sys = DK_sys; W_cm = DK_cm; W_rel = DK_rel

Let's look more closely at the last one:
  W_rel = DK_rel
  W_rel = sum_i int F_i dot dr_ic   [dr_ic = displacement of i relative to cm]
  K_rel = sum_i (1/2)m_i v_ic^2  [v_ic = velocity of i relative to cm]

Careful exmination of the definitions of W_rel and K_rel reveals that
these quantities are simply the total work on the system and total
KE of the system, as measured in the center-of-mass reference frame.

Thus, it seems that W_rel = DK_rel is simply a statement of the W-KE
theorem for the system in the center-of-mass reference frame.

Is it that simple?  Not quite...

This relationship is general (it was proved by subtraction from W_sys=DK_sys
and W_cm=DK_cm), but if the system is not isolated, there may be a net
external force, and thus the center of mass can accelerate.  In such cases,
THE CENTER OF MASS FRAME IS A NON-INERTIAL FRAME.

Does the W-KE theorem apply, as stated, in non-inertial reference frames?

In general, NO.  In proving the W-KE theorem, we used Newton's Second Law:
  dK/dt = ... = ma dot v = F_net dot v = ...
["Proof of Work-KE theorem" in cons-E.txt]

Newton's second law, as normally stated, is only valid in inertial
reference frames.

In this video, we will show
(1) the correction to W-KE theorem for non-inertial frames; and
(2) why this correction is unecessary for the cm frame.

Recall that for a single particle, N2L can be corrected in an
accelerating reference frame as follows:
  F_P,net - m_P a_OG = m_P a_PO  [P=particle; O=observer; G=(inertial)ground]

This was proved by applying N2L F_P,net=m_p a_PG, in an inertial
reference frame, and using relative velocity formulas.  [show page]

The extra force "-m_P a_OG" is added as a correction term to N2L, and
is called an inertial force.  It is not real, but the object behaves
as if it were when examined in an non-inertial reference frame.

If we are observing a system of particles in a non-inertial reference
frame (observer O; a_OG acceleration relative to "ground"), then we
should correct N2L by adding an inertial force
  -m_i a_OG
to each particle in the system.

If we plug this new version of N2L into the proof of W-KE, we find
that
  DK_sys = W_sys + W_inertial
  W_sys = net work done by all REAL forces on the system
  W_inertial = sum int (-m_i a_OG) dot dr_i  [work done by all inertial forces]
W_inertial is the correction term in the W-KE theorem.

  W_inertial = sum int (-m_i a_OG) dot dr_i
             = sum int(ti->tf) (-m_i a_OG) dot v_i dt
             = int(ti->tf) (-a_OG) dot [sum_i m_i v_i] dt
             = int(ti->tf) (-a_OG) dot [M v_cm] dt
             = int(ti->tf) (-M a_OG) dot v_cm dt
Note: dr_i, v_i and v_cm are measured in the non-inertial reference frame

For purposes of work done by an inertial force, we may assume that
the entire combined inertial force for the system [-M a_OG] is applied
entirely at the center-of-mass of the system.

[another way of saying this is that
  W_inertial = W_inertial,cm + W_inertial,rel
==> W_inertial_rel = 0  --  W_inertial is all cm work]

What if our non-inertial reference frame is the center-of-mass frame
itself?

Relative to the center of mass, the center of mass does not
move.  Thus,
  W_inertial = int(ti->tf) (-M a_OG) dot v_cm dt
             = int(ti->tf) (-M a_OG) dot 0 dt
             = 0

No correction for the inertial force is needed in this case, and
  W_rel = DK_rel
can be applied without correction in any situation.

======================================================================

*** System potential energy
**** Conservative internal forces

Previously, we considered conservative vs. non-conservative forces
in the context of a single object subject to a force (presumably
due to an undisclosed second object).
  [see "Conservative forces" in cons-E.txt]

It turns out that the more proper treatment of this concept is to
address it in the context of internal forces within a given system,
which we will now do.


Suppose we have an internal interaction within a system involving one
or more force pairs.
  [diagram showing multiple particles and multiple force pairs]

Examples include:
(1) gravity (pairwise universal gravitation)
(2) Coulomb interaction between electric charges (again, pairwise)
(3) Spring(s) connecting pair(s) of particles
(4) Friction between various objects in the system that are in
contact with one another.

We know that internal forces can do a net non-zero work on the system
(all W_rel, no W_cm).

The question: is this interaction conservative?  What does that mean?


First, some definitions:

(1) The configuration of a system is determined by the position of all
particles in that system.

(2) A process is determined by the (time-dependent) path followed by
the particles of the system as the system transforms from one
configuration to another over a given time interval.

[two-particle system example; with initial and final configuration]

Examples of processes:
* particle 1 moves first (from i->f), then particle 2.
* particle 2 moves first, then particle 1.
* both particles move at the same time.


Consider the work done on all particles of the system by all forces
associated with the given interaction during some process which takes
the system from one configuration to another.

In general, this work will depend on everything:
(1) the initial configuration
(2) the final configuration
(3) the exact process involved

A conservative interaction is one for which the work done by this
interaction depends ONLY on the the initial and final configuration,
and not on the process.

A cyclical process takes the system from one configuration to itself
(initial and final configurations are the same).

An equivalent definition of conservative interaction is one for which
the work done by the interaction (all forces, all particles) is zero
for a cyclical process.

A non-conservative interaction is one for which the work done by the
interaction can depend on the process, or equivalently, there is a
cyclical process for which the work done is not zero.


Examples of conservative interactions are
(1) gravitation (universal, pair-wise)
(2) Coulomb (pair-wise)
(3) Spring(s)

Friction would be an example of a non-conservative interaction.

======================================================================

**** Potential energy

We can associate a potential energy to any conservative interaction
within the system.

The potential energy is defined as a function of the configuration
of the system (i.e., as a function of the position coordinates,
r1, r2, ..., rN for the N-particle system).
  U_sys = U_sys(r1,r2,...,rN)

This potential energy should be regarded as part of the system
(i.e., the total energy of the system includes this potential
energy as a contribution).

The potential energy function can be defined in terms of
potential energy differences, as follows:
  DU_sys = U_sys,f - U_sys,i = -W
  W = work done by interaction forces as the system proceeds i->f

Given that W is process-independent (given that the interaction
is conservative) means that this definition of potential energy
difference is unambiguous.

To calculate U_sys itself, a reference configuration is chosen with an
arbitrary value of U_sys.  Once that choice is made, potential
energies for all other configurations can be determined by the above
relation.

Examples to follow...

======================================================================

**** Example: two objects connected by spring

[see "Two masses connected by spring..." in sys.txt]
[see "Example: spring-mass oscillations with two masses" in sys.txt]

Consider a system with two objects connected by a massless spring
(k) whose natural length is L0.  Assume 1-D motion (x1 < x2).

[diagram]

  Fs21 = -k x21'
  Fs12 = -Fs21 = +k x21'
where
  x21' = x2-x1-L0   distance stretched(+)/compressed(-)

Based on what we know about springs, we expect the spring
interaction to be conservative, and the corresponding
potential energy to be
  U_sys(x1,x2) = (1/2)k(x21')^2 = (1/2)k(x2-x1-L0)^2


Consider the following process where only object 1 moves:
  (x1i,x2) -> (x1f,x2)

Then W2=0 (no motion) and
  W = W1
    = int(x1i->x1f) Fs12(x1,x2) dx1
    = int(x1i->x1f) +k(x2-x1-L0) dx1
    = int(x1=x1i->x1f) +k x21' (-dx21')   substitution
    = -(1/2)k(x2-x1-L0)^2 |_x1i->x1f
    = -U(x1,x2) |_x1i->x1f
    = -[U(x1f,x2) - U(x1i,x2)]


Now consider the following process where only object 2 moves:
  (x1,x2i) -> (x1,x2f)

Then W1=0 (no motion) and
  W = W2
    = int(x2i->x2f) Fs21(x1,x2) dx2
    = int(x2i->x2f) -k(x2-x1-L0) dx2
    = int(x2:x2i->x2f) -k x21' (+dx21')
    = -(1/2)k(x2-x1-L0)^2 |_x2i->x2f
    = -[U(x1,x2f) - U(x1,x2i)]


Putting the two processes together (1 then 2):
  (x1i,x2i) -> (x1f,x2i) -> (x1f,x2f)
  -W = [U(x1f,x2i) - U(x1i,x2i)] + [U(x1f,x2f) - U(x1f,x2i)]
     = U(x1f,x2f) - U(x1i,x2i)

2 and then 1 leads the same result
  (x1i,x2i) -> (x1i,x2f) -> (x1f,x2f)


Now consider an arbitrary process
  (x1i,x2i) -> (x1f,x2f)
where x1(t) and x2(t) are given as functions of time (or
any parameter).

  W = W1 + W2
    = int Fs12 dx1 + int Fs21 dx2
    = int(ti->tf) +k(x2-x1-L0) v1 dt + int(ti->tf -k(x2-x1-L0) v2 dt
    = int(ti->tf) -k (x2-x1-L0) (v2-v1) dt
    = int(ti->tf) -k x21' dx21'/dt dt
    = -(1/2) k x21'(t)^2 |_ti->tf
    = -(1/2) k [x21'f^2 - x21'i^2]
    = -[U(x1f,x2f) - U(x1i,x2i)]


This shows that the spring interaction is indeed conservative, and
the corresponding potential energy is given by
  U(x1,x2) = (1/2)k(x2-x1-L0)^2
as expected.


If the two objects connected by the spring are permitted to move in
3-D, the forces are given by
  Fs21 = -k(|r21|-L0) r21-hat
  Fs12 = -Fs21 = +k(|r21|-L0) r21-hat

This is a central force, so the same techniques can be used as for
the inverse square law (next example).  The result is
  U(r1,r2) = (1/2)k(|r21|-L0)^2

======================================================================

**** Example: two objects interacting via inverse square force

Suppose the two particles in a two-particle system are interacting
with each other via an inverse square force:
  F12 = K/|r12|^2 r12-hat = K/|r12|^3 r12
  F21 = K/|r12|^2 r21-hat = -F12

Is this interaction conservative?

If so, what is the potential energy U(r1,r2) associated with this
interaction?

Imagine that our system moves from one configuration to another:
  (r1i,r2i) -> (r1f,r2f)
via some process given by r1(t) and r2(t).

[see "Work done by internal force..."]

The work done by this internal force pair is given by
  W = W12 + W21
    = int F12 dot dr12
    = int (K/|r12|^2) r12-hat dot dr12
    = -K(1/|r12f| - 1/|r12i|)

To obtain the last step, simply replace r->r12 in the derivation
of the work done by an inverse square force covered earlier.

[see "Work done by a central force..." in cons-E.txt]

Since this work depends only on the initial and final configuration,
and not on the process, the interaction in question is conservative.

Set reference configuration (U=0):
  |r12| = |r1-r2| -> infinity

Then
  U(r1,r2) = U(r1,r2) - U(ref)
           = -W_ref->(r1,r2)
           = +K[1/|r1-r2| - 1/inf]
           = K/|r1-r2|

[see "Potential energy for inverse square forces" in cons-E.txt]

For universal gravity (K = -G m1 m2):
  U(r1,r2) = -G m1 m2/|r1-r2|

For electrical forces (K = q1 q2/(4 pi epsilon_0)):
  U(r1,r2) = (1/4*pi*epsilon_0) q1 q2/|r1-r2|

======================================================================

**** Relate to object view version of conservative/PE

In this video, we will consider the relationship between object-view
and system-view perspectives on conservative forces and potential
energy.

You may have noticed a strking relationship between the object-view
potential energy function (derived previously) and the system-view
potential energy function (derived in the last two videos) for springs
and inverse square forces.  It turns out that this relationship is
quite general.


Suppose we have an object subject to a position-only dependent force
which is governed by a force function f(r), where r represents the
position of the given object (relative to some arbitrary inertial
observer).

If this force is conservative (in the object-view sense), then the
work done by this force as the object moves from ri to rf should not
depend on the path from ri to rf, and so the work can be expressed
in terms of a potential energy function u(r), as follows:
  W = int(ri->rf) f(r) dot dr
    = int(ti->tf) f(r(t)) dot v(t) dt
    = -Du
    = -[u(r(tf)) - u(r(ti))]

The relationship
  int(ti->tf) f(r(t)) dot (dr/dt) dt = -[u(r(tf)) - u(r(ti))]
should be mathematically valid for any path function r(t).


Now suppose that we have an internal interaction between the two particles
of a two-particle system, with the force pair given by:
  F12(r1,r2) = f(r1 - r2)
  F21(r1,r2) = -f(r1 - r2)  [which may or may not equal f(r2 - r1)]
using the same force function f(r).

Q: Is this interaction conservative (in the system-view sense)?

Q: If so, what is the potential energy function?  Can it be related to u(r)?

Well, suppose our system undergoes a process that takes it from
  (r1i,r2i) -> (r1f,r2f)
along some process, specified by path functions r1(t) and r2(t).

The work done by this force pair is given by
  W = W12 + W21
    = int(ti->tf) F12(r1(t),r2(t)) dot v1(t) dt
      + int(ti->tf) F21(r1(t),r2(t)) dot v2(t) dt
    = int(ti->tf) f(r1(t) - r2(t)) dot [v1(t) - v2(t)] dt
    = int(ti->tf) f(r12(t)) dot d(r12)/dt dt
where
  r12(t) = r1(t) - r2(t)
is the relative path of particle 1 relative to particle 2.

Mathematically, that integral can be written in terms of u():
  W = int(ti->tf) f(r12(t)) dot d(r12)/dt dt
    = -[u(r12(tf)) - u(r12(ti))]
    = -[U(r1(tf),r2(tf)) - U(r1(ti),r2(ti))]
where
  U(r1,r2) = u(r1 - r2)

It follows that
(1) the given system interaction
  F12(r1,r2) = f(r1 - r2)
  F21(r1,r2) = -f(r1 - r2)
is conservative, and
(2) the corresponding potential energy function is
  U(r1,r2) = u(r1 - r2)

This can be generalized to an interaction between two particles of a
larger system.  For example, an interaction between particles 2 and
3 of a 4-particle system given by
  F23(r1,r2,r3,r4) = f(r2-r3)   [acts on particle 2]
  F32(r1,r2,r3,r4) = -f(r2-r3)
is conservative, and the potential energy function would be given
by
  U(r1,r2,r3,r4) = u(r2 - r3)

Note that the force functions for the force pair and the potential
energy function do not depend on the positions of the uninvolved
particles (r1 and r4).  Those particles are not involved in the
interaction, and so their positions are irrelevant.

**** Reduction to object-view

Consider again an internal action:
  F12(r1,r2) = f(r1 - r2)
  F21(r1,r2) = -f(r1 - r2)
which is conservative, with a potential energy function given by
  U(r1,r2) = u(r1 - r2)

If r2(t) = 0 (object 2 never moves), then
  F12(r1,r2) = f(r1)
  U(r1,r2) = u(r1)
Furthermore, no work is being done on object 2, since it never
moves.

In such situations, it makes sense to ignore object 2 entirely, and
concentrate the force acting on object 1 and the potential energy
associated with that force, both as a function of object 1's position.

This is precisely how we treated conservative forces and potential
energy in the context of object-view.

For example, with terrestrial gravity, we ignored the earth entirely
and concentrated on the motion of the object itself.

Note that this reduction to object-view only works if the earth
doesn't move, which imposes two limitations:
(1) The earth must remain in place.
(2) We must use the earth reference frame.

If either of these conditions is not met, object-view may fail.

[see "Limitations of object view" in cons-E.txt]


It is interesting to note, however, that system view does NOT fail,
even if the earth is moved (or a different reference frame is used),
since system-view correctly accounts for work done on the earth.

Note that the earth is included in the system when considering
gravitational potential energy in system-view.

Example 1:

Consider the following cylical process:
* Start with earth in place and object at "sea level" (r=R)
   U = -GMm/R
* Now lift object to elevation h
   U = -GMm/(R+h)
   W = -DU = -GMm[1/R - 1/(R+h)]   [work done by gravity on object]
      =approx= -mgh  [if h << R]
* Now move the earth to infinity
   U = -GMm/inf = 0
   W = -DU = -GMm/(R+h)   [work done by gravity on the earth]
* Now move the object back to its original position
   U = -GMm/inf = 0  [earth is still at infinity]
   W = -DU = 0   [no work on object, since F=0]
* Now move earth back to its original position
(object will again be "at see level")
   U = -GMm/R
   W = -DU = +GMm/R  [work done by gravity on the earth]

Overall,
   work done on object = -GMm[1/R - 1/(R+h)]  (1st leg)
   work done on earth  = +GMm[1/R - 1/(R+h)]  (2nd leg + 4th leg)
   overall work done on earth-object system = 0

Conservative force!

Example 2:

Consider the following situation:
  m = 1 kg
  hi = 20m  (release height)
  hf = 0   (hits the ground)
  vi = 0  (measured relative to earth)
  g = 10m/s^2  (makes the numbers easier to deal with)
==>
  dt = 2s  [hi = (1/2)g dt^2]
  vf = -20 m/s  [vf = vi - g dt]

  DU = -mg hi = -200J   [mg = 10N]
  K_obj,i = 0
  K_obj,f = (1/2)m vf^2 = 200J
  Wg = DK = 200J = -DU

Everything checks.

Now consider the same situation, but using a reference frame X
moving down at 5 m/s (v_XE = -5m/s).

Heights will still be measured relative to earth's surface (for
correct U calculation), but kinetic energies will be calculated
from velocities relative to X.

  m = 1 kg
  hi = 20m
  hf = 0
  vi = +5 m/s   [viOX = viOE - vXE]
  g = 10 m/s^2
==>
  dt = 2s  [times don't change]
  vf = -15 m/s

  DU = -mg hi = -200J
  K_obj,i = (1/2)m vi^2 = 12.5J
  K_obj,f = (1/2)m vf^2 = 112.5J
  Wg = DK = 100J   [work done on object]

There is a discrepancy!!

It is resolved by examining the earth during this time.
  v_EX,i = +5 m/s
  M = 6x10^24 kg
  K_E,i = (1/2) M v_EX,i^2 = 7.5x10^25 J   [yikes!!]
  v_EX,f = +5 m/s
  K_E,f = 7.5x10^25 J
  DK_E = 0??

Is that true?  Not quite.  The earth's velocity will actually
change a small amount since there is an upward force on it
(it is attracted to the object).  The change isn't enough to
notice, but it doesn't take much to produce a kinetic energy
change of 100 J.

To calculate DK_E properly, let's look at W_E.

Since v_EX = +5 m/s and gravity is +10N (equal and opposite the
-10N force acting on the 1 kg mass), the work done by gravity on
the earth will be
  DK_E = W_E = F Dx = F v Dt = (10N)(5m/s)(2s) = 100J

Putting it all together:
  DU = -200J
  DK_obj = +100J
  DK_earth = +100J
  DK_sys = +200J

Everything balances!

======================================================================

**** Multiple interactions

Suppose we have an internal interaction which involves multiple force
pairs acting within a given system.

For example, perhaps there are multiple springs connecting pairs of
particles within a given system (rather than just one).

Or perhaps there is an inverse-square interaction between several
pairs of particles.

In fact, universal gravity acts between ALL mass-pairs within a
system, and the Coulomb interaction acts between all non-zero charges
within a system.

Suppose that each force pair (Fij and Fji=-Fij), taken by itself,
is conservative, with potential energy function Uij.

It turns out that the interaction, with all force pairs taken
together, will be conservative, and the corresponding potential
energy function will be given by
  U(r1,...,rN) = sum_ij,pairs Uij(r1,...,rN)

Thus, the potential energy associated with multiple springs would be
given by
  U(r1,...,rN) = sum_ij,pairs (1/2) kij (|ri - rj| - L0ij)^2
where kij is the spring constant and L0ij is the natural length of
the spring connecting i and j.

Likewise, the potential energy for gravity and Coulomb interactions
for a multi-particle system would be given by
  U_g(r1,...,rN) = sum_ij,pairs G mi mj / |ri - rj|
  U_e(r1,...,rN) = sum_ij,pairs [1/4*pi*epsilon_0] qi qj / |ri - rj|

For an N-particle system, there will be N(N-1)/2 such pairs (a
well-known combintorial result.

Proof:

Suppose we have several force pairs, given by force functions
  Fij(r1,...,rN) ; Fji(r1,...,rN)
where
  Fij = -Fji
and that each such force pair is conservative, with a corresponding
potential energy function
  Uij(r1,...,rN)

This means that if the system transforms from an initial
configuration to a final configuration through some process, that
  W^int_ij = Wij + Wji
           = int Fij dot dri + int Fji dot drj
           = int Fij dot drij
           = -DUij
independent of process.

The aggregate work done by the interaction is given by
  W = sum_ij,pairs W^int_ij
    = sum_ij,pairs -DUij
    = -DU
where
  U(r1,...,rN) = sum_ij,pairs Uij(r1,...,rN)

This work done is independent of the process, and corresponds to
the potential energy function above.

If one configuration serves as the reference configuration for every
individual potential energy function Uij (i.e., Uij=0 for that
configuration), then that same configuration will be the reference
configuration for the combined potential energy function.

For example, the reference configuration for the universal gravitation
and Coulomb interactions is the configuration where all particles in
the system are infinitely far apart from each other.  U=0 in that
case.

======================================================================

**** Forces from potential energy function

Given the U(r1,r2,...,rN) function, can we calculate the force on
each particle of system?

It turns out, the answer is yes, in much the same way that we can do
it in the context of object-view.

Consider the situation where particle 1 undergoes a small displacement
(all other particles stay in place):
  r1 -> r1+dr1
  ri -> ri  (i != 1)

The work done by the conservative force pairs associated with this
interaction is given by
  W = -DU
    = -[U(r1+dr1,r2,...,rN) - U(r1,r2,...,rN)]

Because all other particles are at rest, this work is entirely done
by F1 (the net conservative force on particle 1) and can be written
as
  W = F1 dot dr1

It follows that
  F1 dot dr1 = -[U(r1+dr1,r2,...,rN) - U(r1,r2,...,rN)]

If we concentrate on x motion (dr1 = dx1 x-hat), we can proceed as
follows:
  F1x dx1 = -[U(x1+dx1,y1,z1;r2,...,rN) - U(x1,y1,z1;r2,...,rN)]
  F1x = -[U(x1+dx1,y1,z1;r2,...,rN) - U(x1,y1,z1;r2,...,rN)]/dx1
      = -DU/dx1  [all other coordinates fixed]

The same result occurs for y and z:
  F1x = -dU/dx1  [partial derivatives]
  F1y = -dU/dy1
  F1z = -dU/dz1

This can be put together as a gradient:
  F1 = -grad_1 U

"grad_1 U" means the gradient of U, treated as a function of
r1, with r2, ..., rN treated as fixed constants.

Forces on other particles can be calculated in a similar manner,
by considering displacements in each particle respectively.  The
result is
  Fi = -grad_i U


As an example, the gravitational potential energy of a system of
3 point masses is given by
  U(r1,r2,r3) = -G[m1m2/|r1-r2| + m1m3/|r1-r3| + m2m3/|r2-r3|]

The force on particle 1 can be computed by taking the gradient
with respect to r1:
  F1(r1,r2,r3) = -G[m1m2 (r1-r2)/|r1-r2|^3 + m1m3 (r1-r3)/|r1-r3|^3]
which is the net gravitational force on particle 1 due to partciels
2 and 3.

Note that the 2-3 term in U is completely constant in this
calculation, and thus contributes nothing to F1.

Of course, that term contributes to F2 and F3.


The technique of calculating forces from potential energy using
gradients is often used in dynamics simulations of large molecules
(such as proteins and DNA):
* Potential energy functions are developed based on quantum chemistry
principles.
* Forces are computed using gradients, and set to ma for each
particle.
* The time-development of the molecular configuration can then be
worked out by using dynamics simulation.


An interesting consequence of this relationship:

Consider the possibility that all particles in the system undergo
the same displacement, a, and define the following function:
  phi(a) = U(r1+a,r2+a,...,rN+a)

With the help of the chain rule, one can show that
  grad_a phi(a) = sum_i grad_i U(r1+a,r2+a,...,rN+a)
                = sum_i [-Fi(r1+a,r2+a,...,rN+a)]

Since the interaction is completely internal, this sum of forces
will cancel in pairs, resulting in
  grad phi = 0

This means that phi = constant, and so
  U(r1+a,r2+a,...,rN+a) = U(r1,r2,...,rN)

These internal conservative forces do no net work if the system
is moved rigidly as a whole.

This is similar to the result that W_cm is zero for internal forces.

The fact that U is "translation invariant" also implies that it
should generally be possible to express U entirely in terms of
coordinate DIFFERENCES (i.e., r's combine in the form ri-rj).

Consider the gravitational example above.

======================================================================

*** Total mechanical energy

Now that we have correctly defined potential energy, we can define
mechanical for a system:
  E_mech = K_sys + U_sys

K_sys = sum of kinetic energies of individual objects within system

U_sys = sum of potential energies associated with individual
conservative interactions within system

  W_sys = W_ext + W_int
  W_sys = W_ext + W_int,cons + W_int,nc
  DK_sys = W_ext - DU_sys + W_int,nc
-->
  DE_mech = DK_sys + DU_sys
  DE_mech = W_ext + W_int,nc

W_ext includes work done by forces acting on the system from outside
the system, and can increase/decrease the mechanical energy within the
system.
  energy outside of system   <-->   E_mech

W_int,nc includes work done by non-conservative forces within the
system, and can increase/decrease the mechanical energy within the
system.
  other forms of energy within system   <-->   E_mech

Examples include friction, as well as forces resulting from chemical
reactions and other thermodynamic processes (e.g., the car's engine,
the human body, ...).

W_int,cons includes work done by conservative forces within the
system, and does NOT contribute to changes in E_mech.  Changes in
K_sys resulting from these forces are offset by changes in U_sys.
  U_sys   <-->   K_sys   [no net change in E_mech]

Note that if W_ext=0 and W_int,nc=0, then mechanical energy is
conserved:
  DE_mech = 0
  K_sys,i + U_sys,i = K_sys,f + U_sys,f


If our system is the entire universe, then W_ext=0 (there are no
"external forces" in this case, since nothing exists outside this
system).  In that case
  DE_mech,univ = W_nc,univ
and only non-conservative forces can change mechanical energy.

If there are no non-conservative forces, then E_mech is universally
conserved.  We will see how this leads to energy conservation when
we study atomic-view.

======================================================================

**** Two repelling objects released from rest

Suppose we have two objects released from rest and subject to a
conservative and repulsive internal interaction.
  m1 = m2 = 2.0 kg
  U12,i = 40 J
It is assumed that U12->0 as the objects move away from each other.

Also assume that no other (relevant) forces are acting on either
object.

What is v1f and v2f, measured when U12 reaches zero?

----

Examples:

Suppose a massless spring (k=320 N/m) is attached to one of the masses,
and is in contact with the other, and compressed 50cm from its natural
length (x_i=-0.50m).  Then
  U12,i = (1/2)kx_i^2 = (1/2)(320 N/m)(0.50m)^2 = 40 J
U12 will reach zero when the objects move apart and spring expands back
to its natural length.

Suppose the two masses are charged with q1=q2=21.1uC and initially
placed 10cm apart.  Then
  U12,i = (1/4*pi*epsilon_0) q1q2/r12_i
        = (8.99x10^9 Nm^2/C^2) (21.1x10^-6 C)^2/(0.10m)
        = 40 J
U12 will approach zero as r12->infinity.

----

Since the only forces doing work are the internal conservative
interaction forces already noted, we have
  DE_mech = W_ext + W_int,nc = 0 + 0 = 0
and so
  K_sys,i + U_sys,i = K_sys,f + U_sys,f
  K1i + K2i + U12,i = K1f + K2f + U12,f
  0 + 0 + 40J = K1f + K2f + 0
-->
  K1f + K2f = 40J

The 40J of potential energy will be converted to kinetic energy,
and will be shared between the two objects.

Since m1=m2, both objects share KE equally (argue by symmetry,
or using conservation of momentum), we have
  K1f = K2f = 20J
  v1f = v2f = sqrt(2K/m) = sqrt(2*20J/2.0kg) = 4.47 m/s  [4.5 m/s]
Objects will be moving in opposite directions.

Friction and other non-conservativity effects (e.g., spring is not
perfect, EM waves are emitted, spring isn't really massless) will
result in K1f and K2f slightly less than 20J, and the final speeds
slightly less than 4.47 m/s.


What if one object is held fixed, while the other object is free
to move.

The object held fixed must be acted on by an external force, which
holds it in place.

Since F_ext != 0, momentum will not be conserved.

However, W_ext = 0, since the object doesn't move, so mechanical energy
is still conserved.

--> K1f + K2f = 40J

If object 1 is held fixed, then K1f=0 and K2f=40J.  It follows that
  v1f = 0
  v2f = sqrt(2*40/2.0) = 6.32 m/s  [6.3 m/s]

Since m1 cannot move away, m2 will experience its force for a longer
period of time, and will therefore acquire more energy / momentum.


Go back to both objects being released.  What if the masses are not
equal?

Example:
  m1 = 1.0 kg
  m2 = 3.0 kg

Then K1f and K2f may not be equal (symmetry no longer applies).

However, F_ext=0, and so momentum is conserved.

Given the setup, the problem is entirely 1-D.

  m1 v1f + m2 v2f = m1 v1i + m2 v2i = 0
  v1f = -(m2/m1) v2f
  K1f + K2f = U_i
  (1/2) m1 v1f^2 + (1/2) m2 v2f^2 = U_i
  (1/2) [m1 (m2/m1)^2 + m2] v2f^2 = U_i
  (1/2) [m2/m1 + 1] m2 v2f^2 = U_i
  [m2/m1 + 1] K2f = U_i
  [m1/m2 + 1] K1f = U_i 
  K2f = 40J / [3+1] = 10J
  K1f = 40J / [1/3+1] = 30J  [lower mass object receives more KE]
  v1f = sqrt(2*30/1) = 7.75J
  v2f = sqrt(2*10/3) = 2.58J

Note that v1f/v2f = 3/1, consistent with momentum conservation.

Alternative solution:
  v_cm,f = v_cm,i = 0; K_cm,f = 0
  40 J = K_sys,f = K_rel,f = (1/2)mu v21f^2
where
  mu = (1/m1 + 1/m2)^-1 = 0.75 kg
-->
  v21f = sqrt(2*40/0.75) = 10.33 m/s  (>0 based on drawing)
  v1f = v_cm,f - (m2/M) v21f = -(3/4)v21f = -7.75 m/s
  v2f = v_cm,f + (m1/M) v21f = +(1/4)v21f = +2.58 m/s
[see "Two object systems" in sys.txt]


What if there are three objects (or more)?

Then
  K_sys = sum_i K_i  [one term per particle]
  U_sys = sum_ij U_ij  [one term per force-pair]

Assuming that U_sys,i = 40J and U_sys->0 as the objects move
apart, we would find
  K1f + K2f + ... = 40J

If there is adequate symmetry (equal masses, symmetric initial
placement, equal force pairs), then we may be able to conclude
that
  K1f = K2f = ... = KNf
in which case,
  K1f = K2f = ... = 40J / N

If symmetry is not adequate, conservation principles alone may
not be sufficient to solve the problem.  Solution may require
a full Newton's 2nd Law simulation.

For a 3-particle system (moving in a plane), there will be one
energy equation and two momentum equations (one for each component),
but six unknown final velocity components to solve for.

======================================================================

**** Example: two masses + spring (1D)

In this video, we will return to our two-mass-spring problem and
examine it in the context of conservation of mechanical energy.

[see "Two masses connected by spring..." in sys.txt]
[see "Previous problem revisited" in sys.txt]

Situation:

[Force diagrams:
m1: m1g, Fs12, N1G, N1W
m2: m2g, Fs21, N2G]

Two blocks (m1 and m2) are connected by a massless spring (k) whose
natural length is L0 (dist between m1 and m2 centers).  Both blocks
are free to slide (no friction) across a horizontal floor.

At t=0, m1 is pressed against a wall (x1=0) and the spring is
compressed by d (x2=L0-d).  Both are released from rest at t=0.

Findings:

For t <= t0 (while m1 is against the wall):
  x1(t) = 0
  x2(t) = L0 - d*cos(omega*t)
  v2(t) = dx2/dt = d*omega*sin(omega*t)
  omega = sqrt(k/m2)
  t0 = pi/(2*omega) = "period"/4

  xcm(t) = (m2/M) [L0 - d*cos(omega*t)]
  vcm(t) = (m2/M)*d*omega*sin(omega*t)

For t>=t0 (after m1 leaves wall):
  vcm(t) = (m2/M)*d*omega
  xcm(t) = (m2/M)[L0 + d*omega*(t-t0)]

  x21(t) = L0 + A sin(omega'(t-t0))
                ^^^^^^^^^^^^^^^^^^^--- x21'(t)
  v21(t) = A*omega' cos(omega'(t-t0))
  omega' = sqrt(k/mu) = sqrt(k/m1 + k/m2)
  A = d*omega/omega'

  v1(t) = vcm [-cos(omega'(t-t0)) + 1]
  v2(t) = vcm [(m1/m2)cos(omega'(t-t0)) + 1]

Mechanical energy:

The only forces that do any work during any part of this process are
from the spring force pair (Fs12, Fs21), a conservative internal
interaction.

Therefore, mechanical energy should be conserved throughout the motion
of this 2-mass+spring system.

At t=0:
  K_sys = 0  [nothing is moving]
  U_s = (1/2)kd^2  [spring is compressed by d]
  E_mech = K_sys + U_s = (1/2) kd^2

This value of E_mech should remain constant over time.

During 0 <= t <= t0 (m1 is still pushed against the wall):
  K_sys = K1 + K2
        = 0 + (1/2)m2 [d*omega*sin(omega*t)]^2
        = (1/2) m2 d^2 omega^2 sin(omega*t)^2
        = (1/2) m2 d^2 (k/m2) sin(omega*t)^2
        = (1/2) kd^2 sin(omega*t)^2
  U_s = (1/2) k x21'^2
      = (1/2) k [-d*cos(omega*t)]^2
      = (1/2) kd^2 cos(omega*t)^2
  E_mech = K_sys + U_s
         = (1/2) kd^2 [sin^2 + cos^2]
         = (1/2) kd^2

Notation:
  C = cos(omega'(t-t0))
  S = sin(omega'(t-t0))

During t >= t0 (m1 leaves the wall, never to return):
  K_cm = (1/2) M v_cm^2
       = (1/2) M [(m2/M)*d*omega]^2
       = (1/2) M m2^2/M^2 d^2 (k/m2)
       = (1/2) kd^2 (m2/M)
  K_rel = (1/2) mu v21^2
        = (1/2) mu [A*omega'*C]^2
        = (1/2) mu A^2 omega'^2 C^2
        = (1/2) mu d^2 omega^2 C^2
        = (1/2) (m1m2/M) d^2 (k/m2) C^2
        = (1/2) kd^2 (m1/M) C^2
  U_s = (1/2) k x21'^2
      = (1/2) k [A S]^2
      = (1/2) k [d^2 omega^2 / omega'^2] S^2
      = (1/2) kd^2 (k/m2) / (k/mu) S^2
      = (1/2) kd^2 (mu/m2) S^2
      = (1/2) kd^2 (m1/M) S^2
  E_mech = K_sys + U_s
         = K_cm + K_rel + U_s
         = (1/2) kd^2 [(m2/M) + (m1/M)C^2 + (m1/M)S^2]
         = (1/2) kd^2 [(m2/M) + (m1/M)]
         = (1/2) kd^2

Note that once m1 leaves the wall, K_cm can no longer change,
and is no longer "available" for exchange with spring potential
energy.

In such cases, it makes sense to split out K_cm (it is separately
conserved, due to conservation of MOMENTUM), and define
  E_mech,rel = K_rel + U_sys

This represents the mechanical energy of the system as measured in
the center of mass reference frame.

In this case,
  E_mech,rel = K_rel + U_s
             = (1/2) kd^2 (m1/M) [C^2 + S^2]
             = (1/2) kd^2 (m1/M)

The reduced amount of available mechanical energy is responsible for
the decrease in amplitude by a factor of sqrt(m1/M) (from d to A).

[see "Example revisited" above]


Direct calculation of K_sys:
  K1 = (1/2)m1 v1^2
     = (1/2)m1 [vcm (-C + 1)]^2
     = (1/2)vcm^2 (m1 C^2 + m1 - 2Cm1)
  K2 = (1/2)m2 v2^2
     = (1/2)m2 [vcm((m1/m2)C + 1)]^2
     = (1/2)vcm^2 [m2(m1/m2)^2 C^2 + m2 + 2 m2(m1/m2)C]
     = (1/2)vcm^2 [m1(m1/m2) C^2 + m2 + 2Cm1]
  K_sys = K1 + K2
        = (1/2)vcm^2 [m1(1+m1/m2)C^2 + M + 0]   (cross terms cancel)
        = (1/2)vcm^2 [m1(M/m2)C^2 + M]
        = (1/2)M vcm^2 [(m1/m2)C^2 + 1]
        = [(1/2) kd^2 (m2/M)] [(m1/m2)C^2 + 1]
        = (1/2) kd^2 (m1/m2) C^2 + (1/2) kd^2 (m2/M)
        = K_rel + K_cm  (as previously calculated)

[see "Previous problem revisited" in sys.txt; plot v1,v2 vs t]

Note the complex transfer of energy between K1, K2, and U_s.  There
are times when v1 is quite small, and other times when v2 is quite
small.

If we consider the special case where m1=m2, then
  K1=0 at omega'(t-t0) = 0, 2pi, 4pi, ...
  K2=0 at omega'(t-t0) = pi, 3pi, 5pi, ...
  U_s=0 at omega'(t-t0) = 0, pi, 2pi, 3pi, ...

This is "inchworm" behavior.

The energy cycle:
omega'(t-t0)    block 1     block 2       spring
    0           at rest    max speed      relaxed     [picture]
                K1 incr     K2 decr       U_s incr   K2 -> U_s,K1
   pi/2        mid speed   mid speed   max stretched  [picture]
                K1 incr     K2 decr       U_s decr   K2,U_s -> K1
    pi         max speed    at rest       relaxed     [picture]
                K1 decr     K2 incr       U_s incr   K1 -> U_s,K2
  3pi/2        mid speed   mid speed   max compress   [picture]
                K1 decr     K2 incr       U_s decr   K1,U_s -> K2
   2pi          at rest    max speed      relaxed     [picture]
...and repeat...

======================================================================

**** Conservative external forces

When we began studying conservative forces, and, in particular
(terrestrial) gravity, we first adopted an object-view perspective,
and defined the concept of conservative force and potential energy in
terms of the position of the object alone, taking the earth's presence
(and the force it exerts) for granted.

This works well in many cases, but can be problematic if the earth
ever moves.

We then later learned that if we wanted to treat gravity more properly
as a conservative force (using the system-view perspective), we needed
to include the earth in our system and treat gravity as an internal
force.

We then defined the potential energy in terms of the positions of both
the object and the earth, and could now handle situations where the
earth moves.  It is also appropriate to associated the potential
energy with the object-earth system.

However, the earth is rather big place, and we might have reason NOT
to want to include it in our system.  After all, it is likely that
occurances happening elsewhere in the world are of no concern to us.

Of course, we could just go back to object-view and forget about
systems, but this might not work well if there are several objects
that we would like to include in our system, and especially if there
are (other) internal interactions which happen to be conservative.

An example would be to take our 2 block + spring system and place
it on a frictionless incline.


This leads us to the concept of a conservative external force.

To set the stage for this concept, we consider an interaction between
an object (or several objects) within our system and an external
"partner" object which lies outside of the system.  There could also
be several partner objects.

In the case of terrestrial gravity, the partner object would be the
earth itself.

For a given external force to be considered "conservative", the
following must apply:

(1) The interaction between the system and the external partner
object(s) would need to be conservative in the context of a system
which includes the partner object(s).

(2) The partner object(s) cannot move.  This ensures that our system
is effectively moving in a "static force field" generated by the
partner object(s).  This also requires that we use an appropriate
reference frame (one for which the partner object(s) is at rest).

If we were to include the partner object(s) in our system (adopting a
system-view perspective), the potential energy, U, would be defined as
follows:
  U = U(r1,...,rN;r_partner)
  W_tot = W + W_partner = -DU  [independent of process]
  W = work done on system due to this external interaction
  W_partner = work done on partner object

If we assume that the partner object(s) does not move, then we can do
the following:
- delete r_partner from the variable list for U, since it is fixed.
- remove W_partner from the equation, since it is zero.

Thus, U depends only on the system object positions (r1,...,rN), and
  W = work done on system due to this external interaction
    = -DU  (independent of process)


If we have one or more external conservative forces, we may wish
to augment our definition of E_mech to include U_ext (does U_ext
really belong to the system?).
  E_mech = K_sys + U_sys + U_ext
  K_sys = total kinetic energy of system
  U_sys(r1,...,rN) = potential energy associated with internal conservative forces
  U_ext(r1,...,rN) = potential energy associated with external conservative forces

We can also divide W_ext into conservative and non-conservative
parts:
  W_ext = W_ext,cons + W_ext,nc

In this case,
  W_sys = W_ext,cons + W_ext,nc + W_int,cons + W_int,nc
  DK_sys = -DU_ext + W_ext,nc - DU_sys + W_int,nc
-->
  DE_mech = W_ext,nc + W_int,nc = W_nc

External conservative forces would no longer contribute to changes
in E_mech; work done by such forces would be offset by changes in
U_ext.

Suppose we consider our 2 block + spring example on a frictionless
incline.  The forces acting on our system include
(1) Normal forces, which do zero work
(2) Spring forces (internal conservative)
(3) Gravity (from earth) acting on both blocks (external conservative)

There are no non-conservative forces doing non-zero work, and so
mechanical energy, as augmented above
  E_mech = K1 + K2 + (1/2)kx21'^2 + m1 g h1 + m2 g h2
would be conserved during the entire motion.


Some things to note:
  U(r1,...,rN) = U_sys(r1,...,rN) + U_ext(r1,...,rN)
-->
  F_i = -grad_i U  [sum of conservative internal+external forces acting on i]

U_sys(r1,...,rN) is translationally invariant because internal forces
cancel in pairs.

U_ext(r1,...,rN) may not be translationally invariant.

======================================================================

**** The effect of kinetic friction on mechanical energy

Kinetic friction forces, along with other motion-based resistive
forces (such as drag forces), are definitely not conservative.

It follows that the presence of friction within a system of objects
will have an effect on the mechanical energy of the system.

Consider two objects (1 and 2) interacting with each other by a
kinetic friction [diagram].

In this case
  W^fk = W_12 + W_21
       = int fk_12 dot dr_1 + int fk_21 dot dr_2
       = int fk_12 dot dr_12  (as always)
       = int fk_12 dot v_12 dt
       = int |fk_12| |v_12| cos(180deg) dt    (fk_12 always opposite v_12)
       = -int |fk_12| |v_12| dt

Here we are using the fact that fk_12 always points in the (180 deg)
opposite direction as the relative velocity v_12 between the two
objects.

If |fk_12| remains constant throughout the motion (as would be the
case if |N_12| and mu_k12 both remain constant), then it can be
factored out the integral, leading to
  W_fk = -|fk_12| int |v_12|dt
       = -|fk_12| d_12
where d_12 is the total DISTANCE travelled by object 1 in object 2's
reference frame.

Note that this distance can depend on process (not just initial and
final configuration), and will likely not be zero for a closed
process (initial=final), and so kinetic friction is not conservative.

Note also, so long as there is ANY relative motion between the
two objects, that W_fk is always negative.

Since fk is non-conservative, work done by fk contributes to
changes in mechanical energy:
  DE_mech = W_ext + W_int,nc
          = ... - |fk|d
                ^^^^^^^--- contribution to DE_mech from fk
and has a net negative impact on E_mech.

Most of the energy lost from E_mech is converted to thermal energy
within the two interacting objects.  More on this later.

These results also apply to other resistive forces, where f_12
is 180 deg opposite to v_12.


It is important to note that the W_fk we calculated is the
total work done on BOTH interacting objects.  If we focus on
just one of the objects, it is possible that fk can actually
do positive work on that object.

Example: [demo 140b]

In this demo, the book on cart slides forward as cart is moved quickly
(you may recognize this from our discussion of non-inertial reference
frames).

The book starts from rest, and GAINS kinetic energy, due to fk.

However, a more detailed analysis [force diagram] shows that fk
also does work on the cart, and that work is negative.

In fact, the work done by fk on the cart is more negative than the
positive work done by fk on the book, for a net negative total work.
This should be clear based on the fact that the cart is moving faster
than the book (the reason for the drag forward on the book).

As we have proved, this must ALWAYS be the case.

Another example of this effect involves air drag acting on objects.

If the air is moving (i.e., wind is blowing), then objects can actually
gain KE as the air pushes on it and moves it.

However, there is still a net negative work done by air resistance when
the force acting on the AIR (due to the object) is considered.


Calculation example: stacked blocks

Situation [force diagram]:
  m1 stacked on m2  (m2 is a long "runway")
  both slide across a frictionless floor.
  m1 is moving with initial velocity v1i (to the right).
  m2 is at rest initially
  muk = coefficient of kinetic friction between m1 and m2.

Determine the final speed of the m1-m2 combination after their
relative motion ceases, and determine initial and final energies
(various forms) of the system.

  m1 = 1.5 kg
  m2 = 3.5 kg   (M = 5.0 kg)
  v1i = +10 m/s
  v2i = 0
  |fk12| = muk12 m1 g = 2.1 N  (muk12 = 0.143)

Solution #1: Newton's laws and kinematics

Motion:
  a1 = -2.1N / 1.5 kg = -1.4 m/s^2
  v1(t) = v1i + a1 t = 10 - 1.4t
  a2 = +2.1N / 3.5 kg = +0.60 m/s^2
  v2(t) = v2i + a2 t = 0.60t

When does relative motion cease?
  v1(t) = v2(t)
  10 - 1.4t = 0.60t
  10 = 2.0t
  t = 5.0s

Final velocities:
  v1f = v2f = 10 - 1.4(5.0) = 0.60(5.0) = +3.0 m/s

Displacements:
  Dx1 = v1i t + (1/2)a1 t^2
      = 50m - 17.5m
      = +32.5m
  Dx2 = v2i t + (1/2)a2 t^2
      = 0 + 7.5m
      = +7.5m
  Dx12 = Dx1 - Dx2
       = +25m

Initial mechanical energy:
  K1i = (1/2)m1 v1i^2 = 75J
  K2i = (1/2)m2 v2i^2 = 0J
  K_sys,i = K1i + K2i = 75J
  U_sys = 0  (no nonzero-work conservative forces involved)
  E_mech,i = 75J

Final mechanical energy:
  K1f = (1/2)m1 v1f^2 = 6.75J
  K2f = (1/2)m2 v2f^2 = 15.75J
  E_mech,f = K_sys,f = K1f + K2f = 22.5J
  DE_mech = -52.5J

Work done by fk:
  W_fk12 = fk12 Dx1 = (-2.1)(+32.5) = -68.25J
  DK1 = K1f - K1i = 6.75J - 75J = -68.25J  (check)
  W_fk21 = fk21 Dx2 = (+2.1)(+7.5) = +15.75J
  DK2 = K2f - K2i = 15.75J - 0J = +15.75J  (check)
  W_fk = -|fk|d = -(2.1)(25) = -52.5J  (check)

Center of mass / relative motion:
  vcm,i = (m1 v1i + m2 v2i)/M = +3.0 m/s  (coincidence?  no)
  K_cm,i = (1/2)M vcm,i^2 = 22.5J
  K_rel,i = (1/2)m1 v1ci^2 + (1/2)m2 v2ci^2
          = (1/2)(1.5)(+7.0)^2 + (1/2)((3.5)(-3.0)^2
          = 36.75J + 15.75J
          = 52.5J
  K_sys,i = K_cm,i + K_rel,i = 75J  (check)

  vcm,f = +3.0 m/s  (obvious)
  K_cm,f = (1/2)M vcm,f^2 = 22.5J
  K_rel,f = 0J  (no relative motion)

Lost mechanical energy is entirely due to reduction of K_rel from its
initial value to zero.

Solution #2: Conservation principles

Center of mass / relative motion:
  vcm,i = (m1 v1i + m2 v2i)/M = +3.0 m/s
  K_cm,i = (1/2)M vcm,i^2 = 22.5J
  K_sys,i = (1/2)m1 v1i^2 + (1/2)m2 v2i^2
          = 75J
  K_rel,i = K_sys,i - K_cm,i = 52.5J  (or calculate from v_1ci and v_2ci)

No external forces ==> momentum is conserved

  vcm,f = vcm,i = +3.0 m/s
  K_cm,f = 22.5J
  K_rel,f = 0J  (no relative motion remains)
  K_sys,f = 22.5J
  DE_mech = DK_sys = -52.5J
  |Dx12| = d = 52.5J/2.1N = 25m

Individual motions:
  v1f = v2f = vcm,f = +3.0 m/s
  K1f = (1/2)m1 v1f^2 = 6.75J
  K2f = (1/2)m2 v2f^2 = 15.75J
  W_fk12 = DK1 = -68.25J
  W_fk21 = DK2 = +15.75J
  Dx1 = -68.25J / -2.1N = +32.5m
  Dx2 = +15.75J / +2.1N = +7.5m

======================================================================

**** Status of constraint (N/T/fs) and exerted (push/pull) forces revisited

Earlier we attempted to address the issue of constraint forces, such
as normal forces, tension forces, and static friction, and whether
such forces are conservative or not.

[See "Status of constraint..." in cons-E.txt]

The results of our previous discussion were inconclusive.

We will now revisit this issue in the context of system view.


Suppose we have a box sitting on the floor of an elevator.

We already know that if the elevator doesn't move, then the normal
force will do zero work on the box (since it doesn't move either).

However, if the elevator moves upwards, the upward normal force
on the box will do positive work on the box.  However, the
normal force from the box on the elevator (pointing down) will
do negative work on the elevator.

It turns out that whenever we have an interaction force pair
acting on two objects which MOVE TOGETHER, i.e., move in
such way that
  r12 = r1 - r2 = constant
then
  W^int_12 = W_12 + W_21
           = int F12 dot dr1 + int F21 dot dr2
           = int F12 dot dr12
           = 0
since dr12 = dr1 - dr2 = 0 throughout the motion.

The same thing happens with other constraint forces.

Ex: a box lifted by a rope [diagram].

The net work done by tension on the box and the rope together will add
to zero.

Ex: two blocks on top of one other held together by static friction
as one of them (and thus both of them) are pushed to the right [diagram]

The net work done by fs_12 and fs_21 add to zero, even as both blocks
move across the floor.

We also get a net zero work done by an interaction force pair if
the relative motion of the two objects is directed at right angles
to the force.

Ex: a block slides down a frictionless incline, which itself can slide
across horizontal ground.

[see "Block on movable incline plane" in forces.txt]

v12 will be directed down the incline (v1 and v2 will be directed
differently), which is perpendicular to the normal force.  Thus,
  W^int_12 = W_12 + W_21
           = int N12 dot dr12
           = int N12 dot v12 dt
           = 0
W_12 and W_21 separately will not be zero (the incline surely speeds up).

A similar result may occur with tension as long as the rope length
is constant, even if the pivot point is moved.

In general, as long as the constraint is "rigid" (i.e.,
compression/stretch motion can be ignored), the net work done on the
two objects will always add to zero.


What if the compression/stretch motion is not negligible?

Ex: consider two objects undergoing a surface contact collision.

Normally, we would model the normal force as a spring-like
compression force.

If that model holds up, then the net amount of work done by the normal
(spring-like) force on the two interacting objects would depend only
on the initial and final configurations of the two objects (their
initial and final states of compression), and not on the exact process.

[see "Example: two objects connected by spring" above]

The work done would not necessarily be zero, but would be offset by
changes in the potential energies associated with elastic compression
of the two objects.

Note that it may be necessary to treat the objects themselves as
systems, since they are not moving rigidly.

In this context, the normal force would be conservative.

Similar conclusions can be drawn regarding tension forces where
stretching is not negligible, as well as static friction forces where
the two objects may undergo a non-negligible sheer.

[draw the "paralellogram-like" objects interacting via static friction]


As we will explain later, a conservative collision force will
inevitably lead to an elastic collision.

Since inelastic collisions occur, this would suggest that the pure
spring model for normal forces may NOT be adequate in all cases.

Actually, springs themselves can exhibit non-conservative behavior if
they are overstretched, overcompressed, or stretched/compressed too
quickly.

There may be vibrations generated (sound + heat) or permanent
structural damage (bond breaking) which could limit the ability of the
spring to do positive work at the same amount as the negative work
that was done earlier.


The conclusion seems pretty clear: these constraint forces (normal,
tension, static friction) are conservative forces, at least to the
extent that the pure spring model can be applied to them (i.e., no
significant sound/heat/damage generated).

Yet, we struggled coming to any definitive conclusion during our
previous discussion of these forces (object view).

What went wrong?

In the object-view treatment of these constraint forces, the forces
are treated as external forces acting on the object of focus.  The
partner object is not considered.

If the partner object is rigid and does not move (i.e., the floor/hand
doing the pushing, the rope doing the pulling, etc.) then there is no
problem: no work is being done anyway on the constrained object.

However, if the contact surface or rope is moving, then work is being
done on the focus object, and its source is not clear.  There is no
obvious change in potential energy to offset this work.

We can now see that (at least in the rigid case) an equal amount of
opposite-signed work is simultaneously being done on the partner
object (i.e., the surface or rope).

Another issue can come up: suppose the contact point where the
constraint force is applied (surface or rope) does not move, but W_cm
on the object is not zero because the center of mass of the focus
object moves.

In this case, the work done by the constraint force may be zero [see
"External constraint forces" above], but there may be internal forces
(which may or may not be conservative) acting within the focus object
which could change the mechanical energy of that object.

[demo 050.m4v -- I, on a skateboard, push off a wall]

This can give the (erroneous) impression that the constraint forces
are doing non-zero work, and that they may be non-conservative.

In such cases, the constraint force may be redirecting the energy, but
they are not the actual source (or sink) of that energy.

======================================================================

**** Elastic vs inelastic collisions (revisited)

At this point, we will be revisiting elastic vs. inelastic collisions,
this time from the perspective of energy conservation.

[see "The elastic parameter" in cons-p.txt]

Recall that we defined the elastic parameter for a two-body collision,
alpha, by
  |v_12f| = alpha |v_12i|

Collisions can be classified according to this parameter:
  alpha = 0  completely inelastic
  alpha < 1  (partially) inelastic
  alpha = 1  elastic
  alpha > 1  super-elastic / explosive

We also showed that
  K_cm,f = K_cm,i   (this from conservation of momentum: v_cm doesn't change)
  K_rel,f = alpha^2 K_rel,i

Thus, whether a collision is elastic can also be related to changes
(or not) of the total kinetic energy of the two objects (K_sys).

Changes in K_sys can, in turn, be related to the nature of the
collision forces, and whether or not they are conservative.

To simplify the discussion (and remove K_cm from consideration), we
will examine collisions in the center-of-mass reference frame.  In
this reference frame, the two objects will initially be moving
towards each other with zero net momentum.


Consider first a head-on collision (colliding objects move back along
the same line after the collision):

[diagram: before, first contact, full compression, last contact, after
In each case, show objects, their motions, and forces acting]

There will be a certain amount of kinetic energy to start with.

During the first part of the collision, the centers of the two objects
continue moving towards each other as the objects compress.

The force in each case is opposite of this motion, resulting in a net
decrease of the kinetic energy for both objects.

At the midpoint of the collision, both objects will be
(instantaneously) at rest -- zero kinetic energy, and compression will
be maximum.

If the spring-compression model prevails, the decrease in KE of this
system will be offset by increases in elastic PE of the two objects
as they compress.

During the second part of the collision, the two objects will move
apart as the objects uncompress.

The force in each case will be in the same direction as this motion,
resulting in a net increase of the kinetic energy.

Again, if the spring-compression model continues to prevail, the
increase in KE will be offset by decreases in elastic PE of the
two objects as they relax.

If the collision forces are fully conservative, then mechanical
energy will be conserved throughout the collision.
  KE initial --> PE --> KE final

It follows the K_sys,f = K_sys,i, and the collision will be elastic.

However, collisions can be quite violent, and so the perfect
spring-compression model may not prevail.

There may be mechanical energy losses due to friction, sound, or
structural damage to the colliding objects.

In the extreme case (no second part of the collision), the two objects
remain stuck together after the collision, with all of the available
initial kinetic energy converted to other forms of energy (mostly
thermal).

There may also be mechanical energy GAINS if there is a release of
energy (some pent-up PE that is initially present and released along
with the compression PE during the second part; or perhaps a
chemical-based explosion set off by the collision).


Now consider a non-head-on collision in the center of mass frame.

[diagram: before, first contact, full compression, last contact, after
* In each case, show objects, their motions, and forces acting
* Surface normal (n-hat) will be vertical -- J12 assumed parallel to n-hat
* Indicate incident angles, and break velocities into parallel and
perpendicular components (relative to n-hat)]

In this case, it makes sense to separate motion parallel to the
impulse and perpendicular to the impulse, and to use the
coefficient of restitution (e) instead of alpha:
  v12 = v12,par + v12,perp  (vectors)
  K_rel = (1/2)mu v12^2
        = (1/2)mu v12,par^2 + (1/2)mu v12,perp^2    (Pythagoras)
        = K_rel,par + K_rel,perp
  v12,perp,f = v12,perp,i    (no impulse in that direction)
  v12,par,f = -e v12,par,i   (definition of e)
  K_rel,perp,f = K_rel,perp,i
  K_rel,par,f = e^2 K_rel,par,i

At this point, we can employ the same argument.

During the first part of the collision, the centers of the two objects
continue moving towards each other (at an angle to J) as the objects
compress.

The angle between the force and the motion in each case is obtuse,
resulting in negative work done and a decrease in KE for both objects.

At the midpoint of the collision, both objects will still be moving
(because of v12_perp), but v12_par will be zero -- KE will be minimum.

During the second part of the collision, the two objects will move
apart (at an angle) as the objects uncompress.

The angle between the force and the motion in each case is acute,
resulting in positive work done and an increase in KE for both objects.

If the spring-compression model prevails, mechanical energy will be
conserved throughout the collision:
  K_rel,par,i --> PE --> KE_rel,par,f   (e=1)
In this case, the collision is elastic.

If the spring-compression model does not prevail, the collision will
be inelastic, in the same manner as the head-on case.

Note that K_rel,perp (along with K_cm) "goes along for the ride".

It makes sense that the relationship between K_rel,par,f and
K_rel,par,i (which is characterized by e, not alpha) would depend on
the detailed characteristics of the colliding objects, but would NOT
depend much at all on K_rel,perp.


What happens if there is significant friction between the colliding
objects at their physical surface of contact, so that J12 is not
perpendicular to that surface?

We had mentioned before that it may pay to define n-hat in terms of
the surface, and to model the frictional part of the collision force
separately from the normal force part.

In that case, both K_rel,par and K_rel,perp can change, and, in
fact, K_rel,perp will surely decrease due to friction.

======================================================================
